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Abstract 

We give an approach to open quantum systems based on formal deformation quantization. 
It is shown that classical open systems of a certain type can be systematically quantized into 
quantum open systems preserving the complete positivity of the open time evolution. The usual 
example of linearly coupled harmonic oscillators is discussed. 
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1 Introduction 



Attempts at the quantization of open systems, especially dissipative systems, have been made for 
quite some time. Examples can, among many others, be found in [7,11,20]. In particular, some 
approaches to the deformation quantization of genuinely dissipative systems have been conducted, 
see [13,14]. So far, it seems that no successful attempt has been made at a mathematically consistent 
systematic quantization of open systems originating from coupled systems. 

We chose the framework of deformation quantization. The central object of deformation quan- 
tization [3] is the algebra of observables. States are regarded as a derived concept in the sense of 
normalized positive linear functionals on the algebra of observables in the classical as well as in the 
quantum case. The star products used to deform the classical algebra of observables in this process 
are meant to be Hermitian star products. The existence of such star products on the smooth func- 
tions of Poisson manifolds has been proven by [18]. For the special case of symplectic manifolds the 
existence has been proven earlier by [12,16,19], see also the textbook [24] for additional references. 

In the manner of speaking of [6], we get an open system (classical and quantum mechanical) 
by constructing a microscopic model and non-selectively integrating the degrees of freedom of the 
environment. 

As a first step, we give a consistent and general definition of what a classical and quantum 
open Hamiltonian system in the sense of deformation quantization should be relying on the notion 
of completely positive evolutions in both cases. As main result we prove that every classical open 
Hamiltonian system can be deformation quantized preserving complete positivity of the evolution 
map. A by-product of independent interest is the result that for every Hermitian star product 
on a Poisson manifold there is a completely positive map into the undeformed algebra of formal 
series of smooth functions deforming the identity map. Our general formalism is exemplified for 
two coupled harmonic oscillators. 

This article is organized in the following way: In Section [2] a notion of classical open dynamical 
systems in general and the notion of a classical open Hamiltonian system used for deformation 
quantization in particular are defined. In Section [3] we will briefly recall the notions of a Hermitian 
star product and the quantum time evolution with regard to a Hermitian star product. We prove 
in Theorem 13 . 5 1 that for every Hermitian star product one has a completely positive map deforming 
the identity into the formal series of smooth functions with respect to the undeformed product. 
This turns out to be the main tool to show Theorem [3T3J every classical open Hamiltonian system 
can be deformation quantized. In Section 01 as an illustration, we give the standard example of 
the total time evolution of two one-dimensional linearly coupled harmonic oscillators in the setting 
of deformation quantization. Section [5] contains the open time evolutions of a coupled harmonic 
oscillator with respect to states on the bath oscillator corresponding to deformed initial values and 
to KMS states. 

2 Classical Open Dynamical Systems 

There are many ways to specify the notion of open dynamical systems. A fairly general approach 
is obtained as follows: We start with a subsystem whose pure states are described by a smooth 
manifold S and a bath which is described analogously by a smooth manifold B. The combined total 
system has the Cartesian product S x B as space of pure states. 

An open dynamical system is now a time evolution of (pure) states in S x B where we only look 
at the S-part "ignoring" the B-part. More precisely, this is obtained as follows: 

On the total system we specify an ordinary dynamical system, i.e. a vector field X E r°°(T(S x 
B)) with flow : S x B — * S x B. For simplicity, we may assume that the flow ^ is complete, 
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otherwise we have to restrict to certain neighbourhoods in S x B and finite times in the usual way. 
With this assumption, ty t is a smooth one-parameter group of diffeomorphisms of S x B with 

■^-f( = Jo$ t for all teR. (2.1) 
at 

Next we consider the canonical projection maps 

S^SxB^B, (2.2) 

which allow to decompose the tangent bundle T(S x B) into 

T(S x B) = prf TS pr* TB, (2.3) 

where pnf TS and pr B TB denote the pull-backs of the tangent bundles of S and B, respectively. 

Clearly, the map pr s forgets the degrees of freedom of the bath and thus corresponds precisely 
to the idea that we want to ignore the B-part. However, for the time evolution of S we still have to 
specify an initial condition for the bath as well. For the moment, we restrict ourselves to pure states 
and allow for mixed states later on. Thus let x B G B be a point whence we have the embedding 

l Xb : S B x s ^ (x s ,x B ) eSxB, (2.4) 

which is clearly a diffeomorphism onto its image such that pr s o i = id s and pr B o t = x B is the 
constant map. 

Definition 2.1 (Open time evolution, pure case) For any x B G B the open time evolution 

<J>j B : S — ► S of S with respect to the total time evolution * t of S x B and the pure state x B of 
the bath is given by 

cD^p^ov^o^. (2.5) 

Of course, we have to justify this definition and examine some consequences as well as properties 
of <E> t B . First of all, the map 

$ Ib :RxS9(^ s )h^ b (x s )£S (2.6) 

is clearly smooth. However, it does not have the usual properties of an ordinary time evolution: 

For a fixed time t the map <&^ B needs not to be a diffeomorphism, not even for small times. We 
only have the following "evolution property" which easily follows from the one-parameter group 
property of ^>t- 

Proposition 2.2 For the open time evolution we have 

$o B = id s and ^bOM^b)) ^b^j = 
for all x s 6 S, x B £ B, and s, t € R. 

Example 2.3 Let S = R = B and consider the time evolution 

(x s , x B ) = (x s cos(vt) — x B sin(vt),x s sin(z4) + x B cos(^t)) (2-8) 

on S x B. 
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i.) The simplest case is obtained for v £ R being a non-zero constant. Then the open time 
evolution for x B G B is given by 

$t B ( x s) = x s cos(z^t) - x B sm(ut) (2.9) 
which is a diffeomorphism for small t but the constant map for i/t G | + 7rZ. 

nj We can also consider the case where v is a function on S x B depending only on the radius, 
e.g. u(x s ,x B ) = x| + x B . Then ty t is still a one-parameter group of diffeomorphisms and the 
flow lines are still concentric circles around (0,0). However, the points in S x B spin faster 
the further away from (0, 0) they are. Now the open time evolution is 

$t B (x s ) = x s cos((x 2 s + x 2 B )t) - x B sin((x! + s|)t). (2.10) 
E.g. for x B = this gives <&®(x s ) = x s cos(x^t) which yields 

•f(^s)-o (2.H) 

for all t > 0. Since also &® (0) = for all t we see that &® cannot be a diffeomorphism, even 
for arbitrarily small time t > 0. 

From the example we conclude that the open time evolution <I>^ B in general is not a solution to a 
probably time-dependent differential equation on S alone, i.e. in general there is no time-dependent 
vector field X t G r°°(TS) with 

A^B =Xt0 ^ B . (212) 

Nevertheless, this situation of a time-dependent vector field is a particular case of an open time 
evolution as the next example shows: 

Example 2.4 Let Xt G r°°(TS) be a smooth time-dependent vector field on S and let X G 
r°°(T(S x R)) be the corresponding time-independent vector field 



Mr... I) (X t (x s ),^ , j . (2.1:5) 



where we use the splitting (|2.3p of T(S x R) and the canonical constant vector field on the "bath" 
B = R. For simplicity, we assume that X has a complete flow *$>f Then the open time evolution 
for initial condition x B = of the bath is 

$°(x s ) = pr s (tf t (x s ,0)). (2.14) 

But this is precisely the time evolution of the time-dependent vector field Xt, i.e. we have 

±*° t =X t o$l (2.15) 

as an easy and well-known computation shows. Thus the ordinary time evolution of a time- 
dependent vector field can be viewed as a particular case of an open time evolution in the sense of 
Definition EU 
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In view of the yet to be found quantization of open dynamical systems we consider now the 
effect of an open time evolution on the functions C°°(S) as these will play the role of the observables 
later. The following statement is obvious: 

Proposition 2.5 Let x B E B. Then (^ B )* : C°°(S) — ► C°°(S) is a * -homomorphism for every 
t E R and we have 

(**B)* = (id® < y !BB )o* t *opr*. (2.16) 

Here S x : C°°(S) — > C denotes the ^-functional at x B , i.e. the evaluation of a function at the 
point x B . Moreover, \d®S x is the induced map 

id ® ff,. : C°°(S) § C°°(B) = C°°(S x B) — ► C°°(S), (2.17) 

where <g> denotes the completed projective tensor product. Note that the involved Frechet spaces 
are nuclear anyway. 

Though Proposition 12.51 is a trivial reformulation of the definition of it gives a new point of 
view: to this end, recall that a linear functional too : C°°(M) — ► C is called positive if uo(ff) > 
for all functions / E C°°(M). Similarly, we can define a positive functional on matrix- valued 
functions M n (C oc (M)). Having the notion of positive linear functionals we can define positive 
algebra elements by setting that / E C°°(M) is positive if u>o(f) > for all positive linear functionals 
uq. Then it is a true but slightly non-trivial fact that / is positive iff f(p) > for all points p E M. 
The same holds for matrix-valued functions: a function F E M n (C°°(M)) is positive iff F(p) is 
a positive semi-definite matrix for all p E M. Note that in our approach, this is not a definition 
but a consequence of the more algebraic definition. Finally, a linear map <j> : C°°(M) — * C°°(N) 
is called positive if it maps positive functions to positive functions. More important is the notion 
of a completely positive map: eft is called completely positive if all the canonical extensions : 
M n (C°°(M)) — ► M n {C°°(N)) are positive maps for n E N. Clearly, this is the standard definition 
valid for every *-algebra over the complex numbers <D, see e.g. [22] for a detailed exposition and [8, 
App. B] for a discussion of the case of smooth functions. 

Now we come back to our particular situation: while <&* and pr* are canonically given *- 
homomorphisms of the *-algebras of smooth functions and hence completely positive maps, the 
map id (g> S x can also be interpreted as a positive (and in fact completely positive) map which 
coincides with a *-homomorphism t* "by accident". In particular, we can replace the positive 

X B 

functional 5 X ^ by any, not necessarily pure state ojq of C°°(B), that is, a positive linear normalized 
functional loq : C°°(B) — ► C This yields the following, more general definition of an open time 
evolution: 

Definition 2.6 (Open time evolution, general case) For any state loq : C°°(B) — ► <C of the 
bath, the open time evolution of S with respect to the total time evolution ^ t and the state ojq is 
given by 

($£">)* = (id § wo) °#t °prg. (2.18) 

Remark 2.7 Any positive functional ujq : C°°(B) — ► C is actually a positive Borel measure 
with compact support, see e.g. [8, App. B]: for continuous functions this is the famous Riesz 
Representation Theorem, see e.g. [21, Thm. 2.14], which can be shown to extend to the smooth 
setting. Therefore, any state ojq : C°°(B) — ► C is automatically continuous with respect to the 
smooth topology. Thus the map id ®ojq extends to the completed tensor product making the above 
expression in <|2.18j> well-defined. 
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The notation ($^ )* is of course only symbolic as there is clearly no longer an underlying map 
of manifolds. With this definition we shifted the focus to the observable algebra rather than the 
underlying geometry. 

Proposition 2.8 For any state ojq of the bath, the open time evolution ($^ )* : C°°(S) — ► C°°(S) 
is a completely positive map. 

Proof: Since ^1 and pr* are *-homomorphisms we only have to show that id ®ujq is a completely 
positive map from C°°(S x B) to C°°(S). Thus let F G M n (C°°(S x B)) be given and let x s G S. 
Then we have the embedding t : B — ► S x B whence 

S x o (id <g>Lu ) = S x <S> too = ° (S x <8> id) = <^o ° C • ( 2 -19) 

S S S S 

Since i* is a *-homomorphism, the composition ujq o l* is still a positive functional and hence 

x s x s 
a completely positive map. Thus, applied to F*F, we get a positive semi-definite matrix ojq o 

l* (F*F) = 8 o (id <S> too)(F* F). Since this is true for every point x s G S, we have a positive 

% s s 

element (id <§> u )(F*F) G M n (C°°(S)) proving the claim. ■ 

Remark 2.9 Since any positive functional loq : C°°(B) — ► C is actually a positive Borel measure 
with compact support, the map id (§> ujq indeed means to integrate over the bath degrees of freedom 
with respect to a measure specified by ojq. 

Remark 2.10 Note also that in the case of a (5-functional instead of an arbitrary state ujq, the open 
time evolution actually is a *-homomorphism, in contrast to the case of arbitrary states. However, 
in general, (^")* is just a completely positive map without any further nice algebraic features. 

While up to now we have considered arbitrary dynamical systems, we shall pass to more specific 
ones: we assume to have a Hamiltonian dynamics on the total space of the system and the bath. 
In more detail, we choose the rather general setting of Poisson geometry to formulate Hamiltonian 
dynamics. This framework contains in particular any symplectic phase space such as coadjoint 
orbits, cotangent bundles or Kahler manifolds. However, also the dual of a Lie algebra is a (linear) 
Poisson manifold which is important when dealing with symmetries, see e.g. [24, Chap. 3 & Chap. 4] 
for an introduction. 

Thus, let the state space of the system (S,7r s ) and the one of the bath (B,7r B ) be in addition 
Poisson manifolds with Poisson structures 7r s and 7r B . On the total system S x B we choose the 
product Poisson structure 

vr = pr|7r s +pr B vr B . (2.20) 

This means that for functions f s ,g s G C°°(S) and f B ,g B G C°°(B) the factorizing functions / = 
fs ® /b and g = g s ® g B have the Poisson bracket 

{f,9} = {fs, 9s}s ® Ib9b + fs9s ® {/ B , 9b}b- (2.21) 

The dynamics of the total system is given by the Hamiltonian vector field Xh G r°°(T(S x B)) 
with respect to the total Hamiltonian H G C°°(S x B). Recall that the Hamiltonian vector field 
is defined by Xh = {-,H}. In typical situations, the total Hamiltonian contains three parts: we 
have the Hamiltonian H s G C°°(S) of the system alone, the Hamiltonian H B G C°°(B) of the bath 
alone, and an interaction Hamiltonian H I G C°°(S x B) such that the total Hamiltonian is 

H = pr* s H s + pi* B H B + H, . (2.22) 
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Then the total Hamiltonian time evolution is the flow $( : S x B — ► S x B of Xh which we assume 
to be complete for simplicity and analogously to Definition 12.61 the open Hamiltonian time evolution 
with respect to a given state of the bath is defined as follows: 

Definition 2.11 (Classical open Hamiltonian time evolution) The classical open Hamilto- 
nian time evolution of the system S with respect to a total Hamiltonian time evolution &t of S x B 
and a given state loq of the bath is given as the open time evolution 

: C°°(S) — ► C°°(S) (2.23) 

according to Definition 12.61 

Remark 2.12 Again, unless we have special circumstances, the open Hamiltonian time evolution 
is just a completely positive map without any further algebraic features. In particular, there is no 
reason that (<I>^ )* should preserve Poisson brackets, even for ujq = 5 X ^ being a pure state. 

3 Deformation Quantization of Open Hamiltonian Systems 

In this section we will establish the deformation quantized version of the open Hamiltonian time 
evolution. To this end, we recall that a formal star product on a Poisson manifold (M, ir) is an 
associative C[ [ft]] -bilinear multiplication 

oo 

f*g = J2nrc r (f,g) (3.i) 

r=0 

for f,g£ C°°(M)[[ft]] such that Co(f,g) = fg is the undeformed commutative product, C\(f,g) — 
Ci(g, f) = i{f,g} with the Poisson bracket {•,•}, 1 */ = / = /* 1 for the constant function 1, 
and all C r are bidifferential operators [3], see also [24] for a pedagogical introduction. The reason 
that we chose formal star products where a priori no convergence in ft is controlled, is that for this 
situation we have the powerful existence and classification theorems of deformation quantization 
at hand. Physically, of course, one would like to have convergence or at least some asymptotic 
statements. In many examples this is possible but we shall not enter this rather technical issue 
here any further. 

In the sequel, the case where the star product * is Hermitian will be important, i.e. 

7^g = g*J (3.2) 

for all f,g €. C°°(M)[[h]] where ft = ft is treated as a real quantity. This *-involution will be 
necessary to have the honest interpretation of the algebra (C°°(M)[[ft]],*) as observable algebra of 
the quantum system corresponding to the classical system. 

Having the observable algebra, it is natural to define the states in the same way as classically: 
we use positive linear functionals. Now however, we have to specify first what a positive formal 
series should be. Here we can rely on the following definition. A non-zero real formal power series 
a = Yl^Lro ^ a r £ K-[[ft|] is called positive if its lowest non-zero component is positive, a ro > 0. 
In this case we write a > 0. This is a good definition for many reasons: if we view formal 
series as arising from asymptotic expansions then this is what remains from a positive function. 
More algebraically, R[[ft]] becomes an ordered ring by this definition, hence we can rely on the 
rich and well-developed theory of *-algebras over ordered rings, see e.g. [9, 23] for an overview 
and [24, Chap. 7] for an introduction and further references. 
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For star product algebras we can proceed analogously to the classical case and define a C[[/i]]- 
linear functional uj : C°°(M)[[h]] — ► C[[/i]] to be positive if 

w(J*/)>0 (3.3) 

for all / G C°°(M)[[h]]. It can be shown that it suffices to check <$3JS§ for / G C°°(M) without 
higher orders of h. Analogously, we define positive linear functionals for matrix-valued functions 
F G M ri (C 00 (M)[[fi]]) where the star product is extended to matrices in the usual way. Having 
positive functionals we define / G C°°(M)[[ft]] or F G M n (C°° (M)[[h]\) to be a positive algebra 
element if 

u)(f) > and tl(F) > (3.4) 

for all positive functionals u and £1, respectively. Finally, a C [[h]] -linear map 4> '■ C°°{M)\\fi\\ — > 
C co {N)\\fi\\ between two star product algebras on possibly different underlying manifolds is called 
positive if (f> maps positive elements to positive elements. Equivalently, (j) is called positive if 
uj o (p is a positive functional on C°°(M)[[ft]] for all positive functionals uj on C°°(iV)[[ft]]. The 
map 4> is called completely positive if this is also true for arbitrary matrix-valued functions, i.e. if 
4> {n) : M n (C°°(M)[[h]]) — ► M n (C°°(N)[[h]]) is positive for all n G N. Note that even though these 
definitions are in complete analogy to the classical situation, it is nevertheless crucial to have a 
good notion of positive formal power series in ]R[[ft]]. 

Remark 3.1 It is clear that the above concepts generalize immediately to *-algebras A over a ring 
C = R(i) where R is an ordered ring and i is a square root of —1. Even though many of the following 
considerations generalize to this algebraic framework as well, we shall focus on the more particular 
situation of star products. 

Remark 3.2 In the following, completely positive maps will play a crucial role. It is easy to 
see that positive functionals are in fact completely positive maps. Also *-homomorphisms are 
completely positive. Moreover, the composition of completely positive maps as well as convex 
combinations of completely positive maps are again completely positive. Finally, less evident but 
nevertheless true is the fact that the algebraic tensor product of completely positive maps is again 
completely positive. In general, this last statement is wrong for positive maps. 

To describe the positive C[ [ft]] -linear functionals of (C°°(M)[[7i]],*) one first notes that uj is 
necessarily of the form 

oo 

uj = h r u r with linear maps uj t : C°°(M) — ► C (3-5) 

r=0 

Then the positivity uj{f * /) > in the sense of formal power series immediately implies that 
ujo(ff) > classically, i.e. ujq is a positive C-linear functional. This raises the question whether 
every classical state ujq can be "quantized" into a state uj with respect to the star product. In other 
words, we ask whether every classical state is the classical limit of some quantum state. Physically, 
this is absolutely necessary as quantum theory is believed to be the more fundamental description 
of nature. Fortunately, we can rely on the following theorem [10], even for the case of matrices. 
But first we give a definition which shall simplify the further considerations. 

Definition 3.3 (Square preserving map) A C [[h]] -linear map 8 = id + Y^=i h r §r '■ C°°(M)[[fi\] 
— > C°°(M)[[h]} with differential operators S r , S(l) = 1, and §(/) = §(/) is called preserving 
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squares with respect to *, if there are formal series of differential operators D r j : C°° (M)[[h]\ — > 
C°° (M)[[h]\ for r £ No and / running over a finite range (possibly depending on r) such that 

oo 

H7*9) = ^h r ^D^(f)D r j(g) (3.6) 

r=0 I 

for all /,<? G C°°(Af )[[&]]. 

Remark 3.4 It is fairly simple to see that a map preserving squares according to Definition 13.31 
is in fact a completely positive map from the quantized algebra (C°°(M)[[/i]],*) to the classical 
algebra (C°°(M)[[7t]], •) with the undeformed product. 

Theorem 3.5 Given a Hermitian star product *, there exists a globally defined map S preserving 
squares with respect to 

Proof: By [10] we know that for a Hermitian star product ★ on an open subset U C R n there 
exists a map preserving squares with respect to *, denoted by 

oo 

§(7*s) =J2 hr T,^M D ^)- (*) 

r=0 / 

For the Poisson manifold M with star product * we choose a finite atlas. Note that we can always 
find an atlas consisting of dim(M) + 1 not necessarily connected charts. Denote the domains of 
the charts by U a C M. Next we choose a corresponding subordinate finite quadratic partition of 
unity Xa S C°°(M), i.e. supp^o Q U a and J2 a X a Xa = 1- Now let S Q be the locally available maps 
preserving squares with respect to *\u a with corresponding locally defined differential operators 
D r ,i, a - Then we set 

S (/)=Z) a XaXaSa (f\ Ua ) ■ 

Clearly, this gives a globally well-defined formal series of differential operators with §(/) = §(/) 
and 8(1) = 1. Moreover, since the star product is bidifferential, we have (/ * g)\ua = f\u a * 9\u a 
and hence we can apply (*) to obtain 

oo 

HJ*g) = ^^^XaDrJ^ftXaDrj^g). 

r=0 I,a 



Remark 3.6 Recently, a C* -algebraic version of this theorem was obtained for particular strict 
deformation quantizations in [17]. 

The proof of Theorem 13.51 immediately leads to the following consequence. 

Corollary 3.7 For every Hermitian star product -k on a Poisson manifold there exists an equivalent 
star product with the property that every classically positive linear functional ujq is also positive 
with respect to V. 

Proof: This is now easy, as we take a map S preserving squares with respect to *. Then the star 
product f *' g = S(S _1 (/) * §~ 1 (g)) is easily shown to do the job. ■ 
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Remark 3.8 Rephrasing the result from [10] in terms of Theorem 13.51 says that every classical 
positive linear functional loq can be deformed into a positive linear functional with respect to a 
Hermitian star product. Indeed, ujq o 8 will be such a deformation, even universal for all ujq once 
8 is specified. In general, correction terms in higher orders of h are necessary to obtain positivity. 
Moreover, they are by far not unique and neither is the map 8. This is of course to be expected, both 
from a physical and mathematical point of view. Finally, note that each term ujq o § r is continuous 
in the smooth topology, since the classical functional ojq is continuous and the differential operators 
S r are as well. 

After this discussion of states we also need a notion of time evolution for star product algebras. 
Here we can rely on the following facts. For a given Hamiltonian H E C°°(M)[[H]], where we might 
even allow for some /i-dependent correction terms, we consider the Heisenberg equation 

±f(t) = ^[HJ(t)l (3.7) 

for f(t) G C°°(M)[[7j]]. Note that the right-hand side is a well-defined formal power series since 
the commutator vanishes in zeroth order. For simplicity, we assume that the Hamiltonian vector 
field corresponding to the zeroth order Hq of H has a complete flow & f in this case, one can show 
that (|3.7p has a solution for all times with the following properties: There exists a formal series of 
time-dependent differential operators Tt = id + Y^r=\ on M such that 

At = $* t o T t : C°°(M)[[h]] — C°°{M)[[h)} (3.8) 

is a one-parameter group of automorphisms of * with f(t) = Atf being the unique solution of (|3.7j) 
with initial condition /(0) = /. Moreover, At commutes with the commutator [H, and we have 
conservation of energy AtH = H as usual. Finally, if * is a Hermitian star product and H = H a 
real Hamiltonian then At is even a ""-automorphism for each t. For details on this quantized version 
of the classical time evolution we refer to [24, Sect. 6.3.4] and references therein. 

After this preparatory discussion we come back to our original situation of a coupled total 
system S x B. As we already have a nice separation of the total Poisson structure into the Poisson 
structure of the system and the one of the bath, we shall require the same feature also for the 
quantization. Thus, we assume to have Hermitian star products * s on S and * B on B, respectively. 
Then this immediately induces a Hermitian star product * = * s <g> * B on S x B in such a way that 

(C°°(S)p]],* s ) (c~(S x B)[[h]U) (C°°(B)[[ft]],* B ) (3.9) 
are both *-homomorphisms of the involved star products. On factorizing functions we have 

f*9 = (fs*s9s) ® (/b*b0b), (3.10) 

where / = / s ® / B and g = g s g B for f s ,g a G C°°(S)[[h}} and / B , g B G C°°(B)[[h]]. Clearly, 
becomes the first order limit of (13.101) in the commutators. 

Remark 3.9 It will be crucial for our approach that the algebraic structure of the observables 
is a priori given and will stay untouched. The physical interpretation is, that whatever the time 
evolution will be, the way how certain quantities, the observables, are measured is independent 
of any sort of dynamics but a purely kinematical property of the physical system. Thus our star 
products * 3 * s , and * B will be given once and for all and not changed by the open time evolution. 
Note that this is not the only possibility to deal with open systems: in [14] the star product itself 
was modified in order to describe a damped harmonic oscillator. 
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It is now rather obvious what a good definition of a quantized open Hamiltonian time evolution 
in deformation quantization should be: 

Definition 3.10 (Quantized open Hamiltonian time evolution) Let H G C°°(S x B)[[ft]] be 
a Hamiltonian with complete time evolution At and let oj : C°°(B)[[/i]] — > C[[/i]] be a positive 
(D[[?i]]-linear functional. Then the quantized open Hamiltonian time evolution of S with respect to 
oj is 

A? = (id®w) oA t opr* : C°°(8)[[h]] — ► C°°(S)[[h]]. (3.11) 

Remark 3.11 The above completed tensor product is understood order by order in h. Thus we 
have to require that u = Yl^Lo is continuous in each order of h, i.e. each oj t is a continuous 
linear functional with respect to the smooth topology. In view of Theorem 13 . 5 1 and Remark 13.81 this 
seems to be a very reasonable assumption. 

Remark 3.12 Putting Theorem l3.5l Remark 13.81 the existence of Hermitian star products in [18], 
and the existence of the quantum time evolution of Equation (|3.8p together it is easy to see that 
any classical open Hamiltonian time evolution can be quantized into a quantized open Hamiltonian 
time evolution. Conversely, the classical limit of any quantized open Hamiltonian time evolution 
is a classical open Hamiltonian time evolution for the classical limit of the Hamiltonian and with 
respect to the classical limit of the quantum state by construction as Af = (^°)* + 0(h). 

In view of Definition 13.101 it is tempting to believe that the quantized open Hamiltonian time 
evolution Af is completely positive. Indeed, if we would have used the algebraic tensor product in 
(|3.1ip instead of the completed one (g> in every order of H, then this would be a trivial statement: 
the algebraic tensor product of the completely positive maps id and oj is again completely positive, 
and so is the composition with the completely positive *-homomorphisms At and prg. However, 
the crucial point is that the Frechet topology of the smooth functions and the /i-adic topology 
originating from the ring ordering are not very well compatible. In fact, it is not clear whether 
the completed tensor product is completely positive or not. Note that this is rather different from 
the C* -algebraic case where the completed projective tensor product of completely positive maps is 
always completely positive. From that point of view, the following principal result on the quantized 
open Hamiltonian time evolution is non-trivial: 

Theorem 3.13 Let oj be a positive (D[[h]]-linear functional on (C°°(B)[[/i]],* B ) of the form 

oj = uq o 8 (3-12) 

with 8 preserving squares with respect to * B . Then any quantized open Hamiltonian time evolution 
with respect to oj is completely positive. 

Proof: As pr^ and At are *-homomorphisms, the only thing left to show is that id ® oj is completely 
positive. We extend 8 to matrices as usual. For F G M n (C°°(S x B)[[/i]]) we have 

t * B ((id § S)(** * F)) = £V E C B (DrAF)T *s ^ (D r>I (F)) , 

r=0 / 

since the restriction to x B € B commutes with the pointwise products in f)3.6() . Now let jj, : 
C°°(S)[[/i]] — > C[[H]] be a positive C[[ft]]-linear functional with respect to * s . Then for every x B 

t* x (( A1 g§)(F**F))=^rE^fe (D r ,i( F )T*sC (%(f)))eM„(C)i 

B \ B B / 

r=0 I 
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is positive. So if luq is classically positive, we conclude that u>o ° (fi (8> S)(-F* * F) > 0. But 
wo ° (/U (8> §) = \i ° (id ® (wo o §)) = fj, o (id <8> w). Thus, /i o (id (g) ★i 7 ) is positive for all positive 

functionals /i. This implies that (\d®uj)(F* * F) is a positive algebra element for all matrices F 
and hence id <g> lo is a completely positive map as claimed. ■ 

Remark 3.14 The assertion of Theorem 13.131 is actually true for more quantum states than the 
ones of type f)3. 12|) : we will see examples later on in Proposition 15.41 We also note that a possible 
failure of the complete positivity of Af should be seen as an artifact of the rather fine (and not too 
physical) /i-adic topology of formal power series in h. One would expect reasonable behaviour as 
soon as one enters a convergent regime like strict deformation quantization. 

Remark 3.15 In general, the quantized open Hamiltonian time evolution Af is no *-automorphism 
of (C°°(S)[[ft]],* s ). Furthermore, a close look at Equation (|3TTTj) shows that usually A^oAf ^ A% +t 
as expected from a microscopic system. 

Remark 3.16 Using the notions of super manifolds and star products on super symplectic mani- 
folds according to [4, 15] one can easily extend our formalism to this framework. This way, one can 
incorporate spin systems. 



4 Linearly Coupled Harmonic Oscillators I: Generalities 

As an example, consider the well-known linear coupling of two one-dimensional harmonic oscillators. 
We shall describe a one-dimensional harmonic oscillator as a Hamiltonian system (M,tt,H), given 
by M = T*H q ~ R^p, with Hamiltonian H(q,p) = -^p 2 + m %-q 2 , where m, v e IR, + . The Poisson 
bracket is then determined by 

{q,p} = i, {q,q} = o = {p,p}. 

Now let us take S = M = B. The Hamiltonian system (S x B, tt, H) describing the linearly coupled 
identical harmonic oscillators is then given by the smooth manifold S x B ~ R 2 x R 2 , with 

1s>Ps 1b<Pb 

the corresponding Poisson bracket as given by Equation (|2.21|) . In the following, we shall use the 
same symbols q s ,p s , q B ,p B for the coordinate functions on S, B, and S x B, respectively, in order 
to simplify our notation. In the same spirit, we simply write H = H s + H B + H 1 for the total 
Hamiltonian without the explicit use of prjt and pr^ . For the linearly coupled harmonic oscillators 
the interaction term is given by H I = ^(g s — q B ) 2 , with k S R 4 " being the coupling constant. 
Using the new and still global coordinate functions 



9l = +9 B ), Pi = -^(Ps +Pb)> Q2 = -j=(q s -q B ), p 2 = ~^(p s ~ P B 

we can bring the total Hamiltonian to normal form and find the well-known expression 



l /0 , rav 2 r, rav 2 9 2k , 

H= — {p\+pl)+—ql + ^ql with v 2 K = v 2 + —. 4.1 
2m 2 2 m 

The classical time evolution <3?i is known to be a linear map for all t which we can express in matrix 
form as 



2 



/ cos(W) + cosKi) B -^ M - + s -^ 11 cos(vt) - coatvj) _ "M \ 

V / \ ri / rav ' mv K V / V re / rav rav K 

—m{vaia(ut) + v K shx(v K t)) cos(ft) + cos(v K t) — m(vsm{vt) — v K sm(z/ re i)) cos(i4) — cos(^ K i) 

V / V / rav mv K V I \ re 1 rav rav K 

\—m(v sin(z4) — v K mn(v K t)) cos(vt) — cos(v K t) — m{v&m(vt) + v K sm(v K t)) cos(i4) + cos(v K t) J 
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with respect to the global linear coordinates q s , p s , q B , p B . Thus, the open time evolution <I>^ of 
the open subsystem with regard to the state ujq of the bath takes the form 



cos(^) + cos(^) ^ + ^ 



p s J 2 \ — m{vsva.(vt) + v K sin(z/ K i)) cos(z^i) + cos(f K t) 

+ 1 ( uj (q B ) (cos(z4) - cos(z^t)) + wo(p B ) 

2 \-uJo{q B ) m (fsin(i/i) - f K sin(^ K t)) + u; (p B ) (cos(z4) - cos(f K t)) 




(4.2) 



Analogously to the classical case we shall use the normal coordinates in order to simplify the 
computation of the quantum time evolution of the total system. Moreover, it will be advantageous 
to combine the real qi, Pi, q2, and P2 into complex coordinates which will play the role of annihilation 
and creation "operators" later on. We set 

z k = V^Vk + i-^Pk, <lk= 2-4^(z k + Zk), 

_ , . v i and hence Xthtt (4-3) 

z k = ^m^k ~ i^m^Pk p k = ^ (z k - z k ) 

for k = 1,2 and vi = v, V2 = v K . With respect to these global coordinate functions on M the 
total Hamiltonian can be written as H = ^z\Z\ + ^fzi~Z2- For the Poisson brackets one obtains 
{z k ,zi} = = {z k ,zi} and {zk,zi} = j5 k i for all k,l = 1,2. The Hamiltonian for the system will 
take a slightly more complicated form, namely 

2 2 / 

H s = ^ 1 z 1 + ^(z 2 + z 2 ) 2 -^(z 2 -z 2 ) 2 + -^ = (z 1 + z 1 )(z 2 +z 2 )-^^(z 1 -z 1 )(z 2 -z 2 )- (4.4) 

On the other hand, we will also need "factorizing" complex coordinates with respect to the 
original Darboux coordinates on the system S and the bath B. Hence we set 

z s = Jmvqa + i—^p„, z B = Jmvq B + \—l=p 

\ S .^p ' S and _ , y™7 FB ' (4. 5 ) 

z s = ^wq s - i^Ps z B = s/nwq B - i^P B - 

In these coordinates, the Hamiltonians of the system and the bath are given by H s = ^z s z s and 
H B = ^z B z B . The interaction term now reads H Y = (z s + z s — z B — ~z B ) 2 . Again, for the 
Poisson brackets one finds {z k , z{\ = = {z k ,zi} and {z k , ~z{\ = j6 k i for all k, I = s,b. 

After these preparations, we can specify the star product on the total algebra of observables. 
We take the Weyl-Moyal star product on the total system S x B defined by 

f* Wcyl g 

-VV^ iY~ l fi\ l W dr f dr 9 



2 J \ 2 J l Kr-fy h ^u-- -dq^dpi^ ■ ■ ■ dp ir dp h ■ ■ ■ d Pii dq il+1 ■ ■ • % r ' 

(4.6) 

for f,g e C°°(S x B)[[h]], see e.g. [3] and [24, Chap. 5]. 

Remark 4.1 The Weyl-Moyal star product on a flat symplectic phase space R 2n is uniquely 
determined by the requirement of invariance under the affine symplectic group. Under the usual 
quantization map into differential operators it corresponds to the total symmetrization, see e.g. [24, 
Chap. 5] for a detailed discussion. We also note that * Wcyl = ★™ oyl <8> *^ cyl as required by our general 
framework. 
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While the Weyl-Moyal star product is the most natural one with respect to phase space symme- 
tries, it has certain technical disadvantages: when dealing with harmonic oscillators, for technical 
reasons it will be more convenient to employ a Wick star product. Such a Wick star product is no 
longer unique, but depends on the choice of a compatible linear complex structure on the phase 
space which is nothing but the choice of a harmonic oscillator. Therefore, we will have different 
Wick star products adapted to the various harmonic oscillators on hand: either with or without 
the coupling. In detail, one passes from the Weyl-Moyal star product to the Wick star product by 
means of an equivalence transformation explicitly given by 

S = exp(SA) with * = tt£tr k - (") 

k=l 

Then the Wick star product * Wick is defined by /* Wick g = S(S~ 1 (f) * Woyl S' 1 ^)). Explicitly, * Wick 
is given by 

'* 9 -2^— 1^ q z . ... dz . oY- ...dz- { } 

r=0 T - i u ...,i r =l° Z ^ OZ *rOZ ll OZ lr 

for /, g € C°°(E 4 )[[li]]. Alternatively, we ignore the coupling term and use the complex coordinates 
z$, z s for the system and z B , z B for the bath. This gives the two equivalence transformations 

Ss = exp { h dbf) and Sb = exp ( h di^:) ' (4 - 9) 



acting on functions on S and B, respectively. Analogously to (|4.8p we get Wick star products *^ lck 
and *^ lck for the system and the bath, respectively. Since we ignored the coupling terms in the 
definition of the latter two Wick star products, we have 

SytS s ®S B and * Wick + *f ck ® *^ ick . (4.10) 

The total time evolution with respect to * and H can actually be calculated in a much easier way 
than by solving the corresponding evolution equation (|3 . 7j) : we first compute the time evolution 
with respect to the Wick star product * Wlck , which turns out to be simple, and then transform the 
time evolved observables back using S. 

The total time evolution ,Aj Vlck with respect to the Wick star product is determined by 

^r k / = \ [H, AT^fU^ = {AT^f, H} (4.11) 

for / G C°°(S x B)[[h]] due to the fact that H = ^fz\Z\ + ^z 2 Z2- It immediately follows that the 
time evolution is just the classical one, i.e. Af lck = &t ) an d n o higher order correction terms arise. 
But then it is clear that the time evolution with respect to * is given by conjugation with S since 
SH = H + c with a constant c = h u+ ^ K ■ Hence, we have 

A t = S- 1 o o S. (4.12) 

As a consequence we immediately obtain the following result for the open time evolution with 
respect to the Weyl-Moyal star product: 

Proposition 4.2 The deformed time evolution of the open subsystem with respect to the functional 
lu is given by 

A% = (idgwjo^o^o^opr*. (4.13) 
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Remark 4.3 The *-automorphism At obviously restricts to the polynomials Pol (S x B) [h]. Thus, 
being only interested in polynomial observables leads to a convergent formulation of the deformed 
time evolution of the open harmonic oscillator if the quantized state oj used to reduce the total 
dynamics gives a finite order in h for every polynomial on the bath. This will be the case for the 
deformed 5-functionals as well as for the KMS functionals in Section Thence, here we recover 
the usual quantum mechanical formulation including the convergence in h. 

To further illustrate the above situation we compute the open time evolution of some specific 
observables of the system. Here we still allow for a general state lo. 

As a first step we calculate the total quantum time evolutions of the total system for q s and p s . 
To do so, we will have to evaluate the chain of maps (14.131) applied to these observables. First we 
note that 

Sq s =Qs= S ~ 1( ls and S Ps =Ps = S^Ps- ( 4 - 14 ) 
Then the classical time evolution is linear whence applying the transformation S" 1 again does not 
give additional terms. We conclude that 

Atq s =$tq s and A t p s = ®$p s . (4.15) 

For the Hamiltonian H s of the system the calculation is slightly more complicated: First we note 
that applying S yields an additional constant, namely 

SH S =H s + ^(u + u K - —) . (4.16) 
4 V mu K J 

Now the total classical time evolution of H s is quite complicated and can be computed most easily 
from <&* t zi = exp(— itv)z\ and <&* t Z2 = exp(— kV K )z2 and (|4.4p . The remarkable fact is now that 

A$*If s = AH S = \ ( v + u K - (4.17) 
4 V mv K J 

for all t. Thus applying S*" 1 to ^H s gives ^H s minus the same constant as we obtained in (|4.16p . 
We conclude that also for the Weyl star product 

A t H s = $* t H s . (4.18) 

Replacing the complex coordinates and their (simple) time evolution by the original real coordinates 
we get the explicit total classical and hence also quantum time evolutions for q s , p s , and H s 

1 / r n , u / sin(z4) sin(^ K t)\ 

+ \ (cos(, t ) - cosM)), B + - Pb , (4.!9) 

Tfl 1 

A t p s = $tp s = -— (^sin(z^t) + v K sm{v K t)) q s + - (cos(i^i) + cos(^ K t))p s 

TTI 1 

- — (vcos(vt) - v K cos(v K t))q B + - (cos(ut) - cos(u K t))p B , (4.20) 



and 



A t H s = <$>*H S = —(usin{ut) + u K sm(u K t)y + ——(cos(vt) + cos(u K t)y q z s 
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+ — cos (ut) + cos (u K t 2 + — — ^ + —^-L pi 
\ 8m 8 V tnv mv„ ' ' 



i "I* i ./ v ./ \\2 vnv . > \ -> 

+ — (vsm(vt) — v K sm(v K t)) H — (cos(^t) — cos{u K t)) q B 



I 1 / , . . . , 9 mv 2 ( sin(z^i) sin(zA,i) . . 

+ hHcos(^) - cosKt)) 2 + — H \v 



8m ' 8 \ mv mv K 

+ ( — ^(usin(vt) + v K sin(v K t))(cos(vt) + cos(z^t)) 

mv 2 ( sin(^t) sin(zA,i)\ . . . . ,.\ 

+ —r — ^ + (cos(vt) + cos(^t)) q sPs 



4 \ rrw mv K 

(Tfl 
— (vsm(vt) + v K sva.(v K t))(v sin(z^) - v K s'm(v K t)) 

H — (cos(z4) + cos(^ K i))(cos(z^) - cos(z^i)) j q s q B 

( 1. , . , s'm(vt) sinivU) 

+ - - cos (ut) - cos(u K t))( — ^ + v ' 
\ 4 \ mv mv K 

mv 2 , , . . . . /sin(z^t) sin(iA^)\\ 

+ — cos (ut) + cos (i/„t)) — ^ ? s Pb 

4 \ mv mv K J J 

-(cos(i^i) + cos(^ re t))(^ sin(^t) — v K sin(v K t)) 

mv 2 , . . , ,./sin(z4) sin(V K t)\\ 

+ — -(cos i/t - cos (u K t))( + _^Jl2 p g 

+ ( - — (cos(^t) + cos(^ K t))(cos(i/t) — cos(^ K i)) 
\4m 

mv 2 ( sin(^t) sin(z^ K t)\ / sin(vt) sin(z^ K i)\\ 



mv mv K J \ mv mv, 

+ ( - ^(vsin(vt) - v K sm(u K t))(cos(vt) - cos(v K t)) 

mv 2 , , . , ^ / sm(vt) sm(u K t)\\ 

+ — cos (ut) - cos (u K t)) — ^ JLJLi qp 4 .21 

4 \ mv mv K J J 

The reason for transforming the time evolved observables back to the Darboux coordinate functions 
Is i Psi 9b; an d Pb i s n °t j us t an addiction to extensive exercise: It is in these variables where we 
can apply the final map id <g> lo needed for the open time evolution, where a; is a state of the bath 
with respect to *^ oyl . The procedure is very simple: we will have to replace all bath variables by 
their expectation values with respect to oj, i.e. q B is to be replaced by co(q B ), q B p B is replaced by 
u)(p B p B ) et cetera. We will not write down the explicit formulas as these are now obtained from 
flgJID , (|4T20jh and (fOTj) just by copying. 



Remark 4.4 Note that for these observables, the open time evolutions in the classical and quantum 
regime only differ by the (possibly) different expectation values with respect to co and its classical 
limit ujq. In general, we have to expect additional quantum corrections from the total time evolution 
as well. 
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5 Linearly Coupled Harmonic Oscillators II: Examples 



The first example of a state for the bath is a deformation of the 5-functional. Thus, fix a point 
(QboiPbo) m the bath and consider $(q B0 ,p B0 )- For the Weyl-Moyal star product this will no longer be 
a positive functional, see e.g. [24, Sect. 7.1.3]. However, for the Wick star product on the bath 
the 5-functional will be positive without corrections. Thus using the equivalence transformation 
S B we obtain a positive functional 5^ q p ^ j o S B with respect to the Weyl-Moyal star product. 
Note that the equivalence transformation S B is precisely a map preserving squares with respect 
to the Weyl-Moyal star product which is evident from the explicit formula for *^ lck . In fact, this 
was the first example of a map preserving squares which is also heavily used in the proofs in [10]. 
More physically speaking, ^( Qb ^ Pb ) ° S B corresponds to a coherent state localized around the point 

(^bO'^bo)' 

For this particular state we note that for the observables at most linear in q B and p B the 
operator S B does not have a non-trivial effect. Moreover, for the quadratic terms q B , p B , and q B p B 
the operator S B only gives a correction term in first order of H. Explicitly, we obtain 

S B q B = q B , S B p B = p B , S B {q B p B ) = q B p B , (5.1) 

S B ql = ql + h ^^i and S B pl=pl + h^. (5.2) 

From these computations we see that the open time evolutions with respect to 6^ q p ^ ) o S B 
are given by 

8 oS B / 8 \ * 

A 9bo ' Pb ° (q s ) = [^t qB0 ' PB0 ) q 3 , (5-3) 

8 oS B / 8 \ * 

V B0 ' PB ° (p s ) = N>; bo ' Pbo Ps, (5-4) 
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+ — ( - (u sm(ut) - v K sin(z^ K t)) 2 + v \ sax{yt) sin(z^) I + 2v (cos(ut) - cos(^ K t)) 1 . 

16 \ v \ v K 



(5.5) 



Remark 5.1 The classical open time evolutions of q s , p s , and H s in (|5.3p . (|5.4p . and (|5.5p are 
obtained by replacing the functions q B , p B , and their powers in the Equations (|4.19p . (|4.20p . and 
(|4.2ip by their values at q B0 and p BQ . 

Remark 5.2 The deformation of the 5-functional necessary in order to ensure complete positivity 
leads to non-classical components of the open time evolution. 

Next we will study quantized states fulfilling a formal KMS condition, corresponding to "thermal 
equilibrium states" of the bath. 

To this end, we first recall that for every symplectic star product * for C°°(M)[[?i]] there is a 
unique trace functional 

tr:CnM)[[ft]] —►<□[[&]], (5.6) 
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i.e. tr(f-kg) = tr(g*f). Choosing the normalization of tr appropriately one obtains a positive trace, 
see e.g. [24, Sect. 6.3.5] for a detailed discussion and references. For the Weyl-Moyal star product, 
the trace is known to be 

tr(/)=/ f(x)d 2n x, (5.7) 

JR 2 " 

i.e. the integration with respect to the Liouville volume. In fact, it can be shown that in the 
symplectic case the lowest order of tr is necessarily of this form: it is just the integration over the 
whole manifold with respect to the Liouville volume. 

The second ingredient we need is the ^-exponential Exp, as introduced in [3]. Instead of defining 
the exponential function by means of the series, the following approach favoured in [5], see also [24, 
Sect. 6.3.1], will be used. For H G C°°(M)[[h]} one defines Exp(0H) G C°°{M)[[h\] to be the unique 
solution of the differential equation 

^Exp(f3H) = H*Exp(PH) (5.8) 

with initial condition Exp(O) = 1. The classical limit of Exp(i^) is the ordinary exponential exp(i^o). 
The KMS condition for inverse temperature j3 and Hamiltonian H for a C[[ft]] -linear functional as 
formulated in [1,2] in the context of deformation quantization, leads to the following result: up to 
normalization the KMS functional is uniquely determined and explicitly given by 

MKM S (/)=tr(Exp(-/?if)*/) for / € C^(M)[[h}} (5.9) 

see [5] for the proof and [24, Sect. 7.1.4] for more details on KMS functionals. In particular, we 
note that (|5,9p is a positive functional. 

Remark 5.3 Depending on the Hamiltonian H, fi KM s may or may not be normalizable. Whenever 
the Hamiltonian used permits a normalization by rendering the integrations in /% MS (1) well-defined, 
we will denote - — Sts^kms by w K ms and call it a KMS state. 

MKMS {,-*-) 

Before entering the particular example again, we note that in the symplectic case the open 
quantum time evolution with respect to a KMS functional is necessarily completely positive. This 
will follow at once from this proposition: 

Proposition 5.4 Let the system be an arbitrary Poisson manifold and let the bath be symplec- 
tic. Given the KMS functional fJ, KM s with respect to an arbitrary H B G C°°(B)[[h]] and inverse 
temperature j3, the map \6®n KMS : (C^°(S x B)[[7i]],*) — > (C^(S)[[h]], * s ) is completely positive. 

Proof: We choose a map 8 for the bath whose existence is guaranteed by Theorem l3.5l In the proof 
of Theorem l3.13l we have seen that for every positive C[[fr[] -linear functional fi : M n (Cg°(S)[[?i]]) — > 
C[[ft]] the combined map 

//§ § : (M n (CHS x B)[[ft]]),*) — » (C °°(B)[[ft]], ■) 

is positive. It follows that for F G M n (Co°(S x B)[[ft]]) the function (fj,® S)(F**F) is at every point 
x B £ B either a formal series with positive lowest order term or zero. To avoid trivialities, assume 
that (^<8> §)(F* -kF) is not identically zero. Let ro be the minimal exponent with Qu<g) §)(F* *F) = 
h r °a ro + • • • and a ro > not identically zero. By continuity, there is an open subset [/CB with 
a ro (x B ) > for x B G U. But this implies that (/%ms o§ _1 )o(/i® §)(F* * F) = h r °b ro + ■ ■ ■ with 
b ro > since the zeroth order of S is the identity and the zeroth order of /%ms is the integration 
over all of B. Since [i is arbitrary and using 

(/Jkms o § _1 ) o (/i g 8) = // o (id § /Jrms), 
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this shows that (id ® ^kms)(F* * F) is a, positive algebra element in M n (Cg°(S)[[^]]) with respect 
to * s . ■ 

Back to our specific example, we consider the harmonic oscillator as the Hamiltonian H B E 
C°°(]R 2 ) and the Weyl-Moyal star product * B as before. In this case, the star exponential of H B 
has been computed explicitly by [3]. One has 

Exp( -^ ) = exp (-^ tmh {*¥)) (510) 

for (5 > and v > 0, which is a well-defined formal power series in H. Note that in [3] the 
exponential Exp(j-H) requires a convergent setting due to the h in the denominator. In our case, 
the situation is much simpler. In fact, differentiating (|5,10p with respect to (3 gives the defining 
differential equation (|5.8p right away. 

As in the textbooks on statistical mechanics, we can now calculate the partition function Z as 
the normalization factor of the KMS state on the bath by formally calculating Gaussian integrals. 

Proposition 5.5 The normalization factor fi KMS (l) is explicitly given by 

exp [_Mi£ 

ta(1) = M r^) e «' (5 - n) 

The partition function is the formal Laurent series 

( hfiu 

GXp I — 

z = 1 (\'r\ e ( 5 - 12 ) 

1 — exp \ — ripv) 

The crucial point is that /x KMS (l) has a well-defined classical limit while Z has a simple pole at 
h = 0. Therefore, we can use this normalization factor to obtain the well-defined KMS state 



>(/) = ^ / E M-PH B )* B fdq B dp B (5.13) 



for / E C°°(B)[[^]] such that the integral (|5. 13j) is convergent order by order in h. Note that the 
inverse of 2nhZ is again a well-defined formal power series. 

As for the ^-functional, we shall now compute the open quantum time evolution of the observ- 
ables q s , p s , and H s also with respect to the KMS state co K ms- To this end, we need the expectation 
values of q B , p B , q B , p B , and q B p B in order to evaluate (I4.19p . (I4.20p . and (I4.2ip . 

Lemma 5.6 One has the following expectation values 

Wkms(<?b) = ^kms(Pb) = ^kms(^bPb) = 0) ( 5 - 14 ) 

Wkms(?b) = x , and Wkms(Pb) = 777 \ , ( 5 - 15 ) 

2mi/tanh(M^) 2tanh(^ 



which are formal power series in <C[[h]]. 

Proof: This is of course textbook knowledge. Nevertheless, we sketch the computation in order 
to illustrate the star product formalism used. The first observation is that the trace functional tr 
for the Weyl-Moyal star product has the remarkable feature 

tr(f*g)=tr(fg), 
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see e.g. [24, Ex. 6.3.33]. Strictly speaking, one of the functions has to have compact support. 
However, if one is the Gaussian Exp(— (3H B ) then the rapid decay allows to perform the integrations 
by parts also for observables like polynomials. Thus we can use this feature to simplify co KMS (f) 
considerably for the above observables. Since Exp(— (3H B ) is just a Gaussian we are left with the 
well-known computation of some Gaussian integrals. ■ 

Using these expectation values, we can apply the general formulas (I4,19j) . (|4.20j) . and (|4.21j) and 
substitute there the observables q B , p B , q B p B , q 2 , and p 2 by their expectation values with respect 
to cl>kms • This then gives the open time evolutions of q s , p s and H s . Remarkably, many terms 
disappear thanks to the vanishing of (I5.14p . In detail, we have 

, „„„ 1 , . , . If sm.(vt) sin(f K i)\ 

^ KMS 9 S = » (cos(^) + cos(^t)) q s + - — ^ + p s , (5.16) 

Tfl \ 

^ KMS Ps = -y (vsm(vt) + v K sm{y K t)) q s + - (cos(z/i) + cos(v K t)) p s , (5.17) 

(Tfh TTVlJ^ \ 

— (i/sin(i/i) + v K sm(v K t)) 2 + — -(cos(i/i) + cos(v K t)) 2 j q 2 

(1 . . , , u , mv 2 fsvaivt) sm.(v K t)\ 2 \ 9 

— cos (i/t) + cos (u K t 2 + — - — ^ + V > pi 
8m 8 \ mv mv K J I 

+ ( — —(vsin(yt) + v K sm(v K t))(cos(vt) + cos(^ K t)) 

mv 2 f s'm(vt) sin(z^)\, , . . ,A 

+ ~r — ^ + cos ^ + cos (v K t)) q sPs 

4 \ mv mv K J J 

— {ys\u.{vt) — v K sm.(v fi t)) 2 + 2^(cos(z^i) — cos(z^t)) 2 

v \ <$K 

+ v(sm(vt) - — sin(^)) 2 I (5.18) 

"« /I6tanh(^ 
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